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Single-electron computing without dissipation 
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A possibility to perform single-electron computing without dissipation in the array of tunnel- coupled quantum dots 
is studied theoretically, taking the spin gate NOT (inverter) as an example. It is shown that the logical operation can 
be realized at the stage of unitary evolution of electron subsystem, though complete switching of the inverter cannot 
be achieved in a reasonable time at realistic values of model parameters. An optimal input magnetic field is found as 
a function of inter-dot tunneling energy and infra-dot Coulomb repulsion energy. 



Recent advances in the fabrication of nanometer scale quantum dots open an opportunity of practical implemen- 
tation of the idea to use the states of a quantum system for data coding and processing Q. For example, the spins 
of individual electrons can be viewed as the bits of information: logical one (zero) corresponds to "up" ("down") 
direction of electron spin at a given quantum dot. Spin gates (elementary sets of quantum dots performing particular 
logical functions) have been discussed by Bandyopadhyay et al || and later investigated theoretically by Molotkov 
and Nazin || and by Krasheninnikov and Openov [Q. If occupation/unoccupation of a quantum dot by a single 
electron is viewed as a bit 1/0, one has charge gates, various kinds of which have been studied, e.g., by Lent et al 
Q and by Nomoto et al In arrays of quantum dots, the quantum tunneling of electrons between adjacent dots 
and/or Coulomb interaction of electrons with each other play a role of "wiring", resulting in the signal propagation 
T^j- from dot to dot. 

(f) ■ The operation of spin and charge gates is based on the principle of "ground state computation" According 
to this principle, upon the influence of external source on input dots of a particular gate, the electron subsystem 
changes to a new ground state. The final spin or charge configuration reflects the result of "calculation". This result 
can be read from output dots of the gate. The quantum dot gates are believed to possess high-speed performance as a 
consequence of extremely fast switching between different electron ground states. However, the switching rate, being 
dictated by dissipation processes, is not known a priori. For a gate consisting of a few quantum dots and operating at 
sufficiently low temperatures, the switching rate may appear to be rather small (10 6 to 10 9 s _1 ||), thus slowing down 
the computation. Hence, in studies of the potential of quantum dot arrays for high-speed single-electron computing 
£H . one should give special attention to inelastic relaxation processes. 

An alternative way has been discussed recently by Bandyopadhyay and Roychowdhury [jr). They explored the 
dynamic behavior of the simplest spin gate NOT (inverter) and found that there exists optimal input signal energy to 
achieve its complete switching in the absence of inelastic relaxation. However, it remains unclear how adequately the 
results obtained in depict reality since the authors of Ref. Q used the Heisenberg model to describe the correlated 
electrons in quantum dots. Meanwhile, it is well known that this model is just a limiting case of more realistic 
Hubbard model || and cannot be used if the inter-dot electron tunneling energy V is of the order or greater than the 
intra-dot electron repulsion energy U . It is instructive to study the broad range of V and U values in order to see if 
the conclusions of Q] reflect the basic physics or are just a consequence of using a particular theoretical model. 

In this paper we study the unitary evolution of electron subsystem in the spin gate NOT (inverter) making use of 
the Hubbard model with arbitrary values of V and U. The inverter consists of two closely spaced quantum dots (A 
and B) occupied by two electrons || , ||] , JjJ. One of two dots (say, the dot A) serves for writing the input signal to 
the gate by the action of the local magnetic field Ha- The second dot (B) is the output. At Ha = the ground state 
of the inverter is the entangled state with zero magnetic moments at both dots A and B. The logical function NOT 
is realized if at Ha ^ magnetizations (i.e., spin projections) of dots A and B have opposite directions. Complete 
switching of the inverter is said to take place if spin projections are saturated (S z a — 1/2, S z b = —1/2 or S z b = 1/2, 
S z a — —1/2, where S Z i — (S Z i) — (n^ — rii|)/2, hi being the operators of particles number at dots i = A, B). Upon 
complete switching, magnetizations of both dots reach the maximum absolute value gfiB, where g is the Lande factor, 
Hb is the Bohr magneton. We stress that the ground state at any finite value of Ha is organized in such a way that 
\S z a\ < 1/2 and \S z b\ < 1/2 |J, 0). Hence, the complete switching of the inverter cannot be achieved through its 
inelastic relaxation to a new ground state. 

The Hubbard Hamiltonian for the inverter has the form 

H = -V^2(a\ a a Ba + a,g a a A<J ) + Uh A ^fiAi + Un B ^h B i - gnBH A ^2h A aSign(a), (1) 
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where the quantities V, U, Ha, g, and [Ib are defined above, other notations being standard (see, e.g., 0). Here we 
assume that each dot has one size-quantized level with on-site potential Eq = (i.e., all energies are measured from 
£o). In what follows, we set gfiB = 1. 

The complete orthonormal set of inverter eigenstates is formed by two-electron basis states |1) = | T>l)i |2)H l,T), 
|3) = | 11, 0), |4)=|0, 11), |5) = | t) T)j |6) = | |, |), where, e.g., the notation | f, I) denotes the state with up-spin electron 
at the dot A and down-spin electron at the dot B, the notation | Tl>0) denotes the state with two (up-spin and 
down-spin) electrons at the dot A and no electrons at the dot B, etc. The magnetic moment of the electron with 
up-spin polarization is oriented along the direction of the local applied magnetic field Ha- 

At Ha — 0, the ground state eigenvector of the Hamiltonian (1) is 




16^2 _ u JlP + lW 2 -U \ 

13) + V \ v 14} . (2) 



The corresponding eigenenergy is Eq = (U — \/U 2 + 16T^ 2 )/2. In the ground state we have (\f r o|5' Zj 4|^ r o) = 
(^oPzbI^o) = 0. We suppose that at t < the system is in its ground state. 
If the local external magnetic field is applied at time t = 0, then the wave function at t > is 

6 

*(*) = exp(-iE k t/h), (3) 

k=l 

where ^ k and Ek (k = 1 — 6) are eigenvectors and eigenenergies of the stationary Schrodinger equation 

HV k = E k V k . (4) 
The coefficients A k should be found from the initial condition ^(t = 0) = ^5>q. It is convenient to write ^f k as 
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$ k = ^B kn \n). (5) 

n=l 

Then 

*(«) = £/n(t)|n>, (6) 

71=1 

where 

6 

f n (t) = A k B kn exp(-iE k t/h). (7) 
fe=i 

The probability to find the system in the basis state \n) at time t is p n (t) = \f n (t)\ 2 - 

At arbitrary values of V, U , and Ha the eigenvalue equation (4) reduces to the algebraic equation of the third 
power in E k . The resulting analytic expressions are too cumbersome for analysis, so it is more convenient to solve 
the equation (4) numerically. Before proceeding to the results of these calculations, let us consider the limiting case 
U = which physically corresponds to U « V (closely-spaced large-sized dots |J). 

At U = we have rather simple equations for the probabilities p n (t): 

Pl (t) = -(l+ \ HAV 9 sin 2 M/2)V , p 2 (t) = -(l- \ HaV 9 sin 2 M/2) 
fi\ ! 4 I ^ H 2 + W 2 V 1 ' J ' F K ' 4 I H 2 +W 2 y ' ' 



1 / 16H 2 A V 2 



sin 4 (w*/2)J , p 5 (t)=pe(*)=0, 



Ps(t) = Pi(t) = l[ 1 ~ {H 2 + W 2 f sin 4 (^V2) ) , P5 (t) = Pe (t) = 0, (8) 



where w = \J H 2 A + 4V 2 /U. From (6) and (8) it is straightforward to find that 

S zA (t) = (*(t)|& A |*(t)) = -S zB (t) = (pi(t)-p 2 (t))/2 = H l H +l V 2 sin 2 M/2). (9) 
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From (9) we see that the spins S z a and S z b are oppositely directed at any time t according to the physical truth 
table of the logical gate NOT ||. For the sake of definiteness, let us consider the case Ha > 0. In this case S z a 
is always positive and peaks at to — tt/lu. Moreover, a complete switching, S z A{to) = 1/2 and S Z B(to) — —1/2, is 
achieved at Ha/V — 2. The dependence of S z A(to) on Ha/V is shown in Fig.l. We stress that this dependence 
has been obtained by us in the weak coupling limit of the Hubbard model. Nevertheless, it is analogous to those 
calculated in |Q within the Heisenberg model (i.e., in the strong coupling limit of the Hubbard model, U » V) with 
an exception that in the Heisenberg model lo = U H A + 4J 2 /?i, to — tt/2uj, and S z A(to) reaches a maximum value of 
1/2 at Ha = 2 J |?J, where J is the antiferromagnetic exchange energy (for two-site cluster J = V 2 /U at U » V). 
Hence, one may expect that complete switching of the inverter can occur at an arbitrary ratio of U to V . 

To check this hypothesis, we calculated numerically the dependencies of S z a on t and S Z A(to) on Ha/V at different 
values of U/V, where to is generally defined as a time of the first maximum at the curve S z A{t), to being a function 
of Ha/V and U/V. The curves of S z a{Pq) versus Ha/V are shown in Fig.l for several values of U/V. One can 
see that increase in U/V first results in decreased height of the maximum at the S z A(to) versus Ha/V curve. At 
U/V > 2 the height of this maximum increases again, but doesn't reach the saturated value 1/2 at finite U /V, though 
S z A(to) — * 1/2 if U/V — ► oo (this corresponds to the Heisenberg model and agrees with the results obtained in (t]]). 

It seems that complete switching of the inverter cannot be achieved at realistic values of U/V ratio, i.e., at U/V ^ 
and U/V ^ oo. Note, however, that at arbitrary values of Ha/V and U/V the function S z A(t) is not periodic in 
time since it includes several harmonics with different frequencies and amplitudes. Hence, in principle, the value of 
S z a — 1/2 can be achieved at some longer time. But this case is of no interest for us since we should like not only to 
reach the maximum permissible value of S z a, but to do it in as short as possible switching time. 

However, an impossibility to achieve the complete switching of the inverter doesn't imply the impossibility to 
perform the logical operation NOT at the stage of the unitary evolution. One should just to "read" the signal at a 
time when S z a(b) has a large absolute value, e.g., at a time to- Indeed, S z A(t ) > 0.45 at any value of U/V (see 
Fig.l). Hence, the error probability P err = 1 — pi(t ) (i.e., the probability to read the "wrong" signal S z a = —1/2 or 
S z a — at a time to) is less than 0.1. Our calculations showed that at U /V « 1 and at 'optimal' (for a given U/V) 
value of Ha/V the time to is of the order of h/V, i.e., to ~ 10 _13 s for V w 10 meV. The value of to increases as U 
increases and reaches ~ 6U/V at U/V = 10. The limiting value of to at U >> V is to — nhU / ^s/TV" 1 , in accordance 
with |Q. Thus, to speed up the calculation, we should have small U and large V . If the shape of a single quantum 
dot is a cube with side length a and the distance between the quantum dots is d, then V decreases exponentially 
in both d and a, while U is roughly inverse proportional to a and is almost independent on d 0, Hence, small 
values of d and a favor short switching times to (the values of U and V can be calculated numerically for a given set 
of geometrical parameters of quantum dots array and for a particular choice of semiconducting materials 0). 

On the other hand, the 'optimal' value of Ha increases with V as H° A pt = 2V at U = (9) and as H°f = 2V 2 /U at 

U/V » 1 (see also 0). The product H A pt to is of the order of h at any U/V. Hence, one can have a realistic value of 
H opt < 1 Tegla Qnly at the 

expense of increasing to up to w 10 11 s. Nevertheless, this value of to still remains several 
orders of magnitude smaller than characteristic times of inelastic relaxation . 

In summary, the switching of the spin gate NOT (inverter) at the stage of unitary evolution is much faster than 
through relaxation to a new ground state. The switching time can be reduced down to 10 _11 s through proper choice 
of quantum dots geometrical parameters and local external magnetic field, with the error probability less than 0.1. 
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Figure caption 

Fig.l. The maximum value of S z A(to) versus Ha/V at different U/V. 
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